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This paper describes a novel fuzzy tracking control for permanent magnet syn- 
chronous motors (PMSM) using an extended Kalman filter (EKF) and distur- 
bance observer (DO). The goal is to create a robust controller able to drive the 
system’s states to track a virtual reference model and provide a low disturbance 
effect on the synchronous machine. First, the PMSM is represented using a 
fuzzy model Takagi-Sugeno (T-S) attended by the load torque variation. Next, 
To simplify the construction of a virtual reference model and nonlinear track- 
ing control, a fuzzy tracking control based on virtual desired variables (VDVs) 
is proposed. Using this concept, a two-stage design procedure is developed: 
i) determine the VDVs using the desired output and the load torque, which can 
be estimated using DO and EKF and ii) calculate the fuzzy controller gains by 
solving linear matrix inequalities (LMIs). The efficiency of the suggested strat- 
egy is eventually shown through simulation results. 
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1. INTRODUCTION 

During the last years, the high efficiency, power/weight and torque/inertia ratios of permanent mag- 
net synchronous motors (PMSM) have attracted increasing interest in industrial drive applications. Due to the 
inherent nonlinearities and external disturbances, their analysis and control is a difficult task. Thus, the lin- 
ear control approach cannot guarantee satisfactory performances. To avoid the difficulties connected with the 
design of the PMSM controller, many proposals have been developed in the previous three decades, e.g. Di- 
rect torque control, adaptive control, neural network control, feedback linearization control, and sliding mode 
control have been considered respectively in [I]-[5]. 

Recently, many schemes based on the Takagi-Sugeno (T-S) fuzzy model have been proposed [6]. 
Using a T-S fuzzy-model-based, a complex dynamic system can be described using an averaged sum of a set 
of local linear subsystems, thus providing a systematic framework for modelling, analyzing, and controlling 
nonlinear systems. Using the Lyapunov approach, stability analysis is carried out, where the control issue is 
formulated into linear matrix inequalities (LMIs) problem [7]. 

The problem of tracking control of PMSM using the T-S fuzzy model and a novel concept namely 
virtual desired variables (VDVs) has been studied by many authors [8]-[[11]. Based on this notion, the challenge 
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of tracking control can be transformed into a simple stabilization problem which has been discussed in many 
works using the know parallel distributed compensation (PDC) technique |12]. Moreover, the control law and 
the reference model can be efficiently designed using generalized kinematics [13]. 

However, the control of PMSM using T-S fuzzy models and the VDVs concept has been discussed 
in many works without considering the influence of external disturbances. In practical industrial applications, 
the PMSM is usually affected by disturbances; their presence diminishes the tracking control performance. 
Hence, many works have dealt with the T-S fuzzy control of nonlinear systems with disturbances using several 
techniques, especially those based on H-infinity tracking control [14], [15]. Thus, the first idea is if possible to 
combine the concept of VDVs with H-infinity performance. In this case, it can design a fuzzy controller that 
allows benefiting from the first one, the facilitated design of the reference model and the control laws, and from 
the second one, minimizing the effect of the disturbances on the system. However, this method is not usually 
effective when disturbances are essential, as in the case of the PMSM system. 

This work develops an advanced robust fuzzy tracking scheme for a PMSM. To this end, the PMSM 
is represented by the T-S fuzzy model. In order to simplify the design of nonlinear control law and a reference 
model, a fuzzy tracking control based on the set of VDVs is then constructed. Instead of minimizing the effect 
of the disturbance on the system, the disturbance becomes a part of the reference model to compensate entirely 
for its effect on the motor. However, the problem here is that the disturbances are usually in an unknown 
state, or its measure is a difficult task as in the case of considering motor when the load torque is considered 
an external disturbance. To solve this problem, the load torque is estimated using a disturbance observer 
and extended Kalman filter, as outlined in Figure 1. The stability of the augmented system is analyzed by 
Lyapunov’s approach, which may be expressed as LMIs problem. An algorithm to calculate the controller 
gains according to the desired system characteristics is presented and discussed. Matlab/Simulink simulations 
of the acquired findings demonstrate the great performance of the suggested control approach. 
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u 
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Figure 1. Proposed fuzzy tracking control scheme 


2. PMSM MODEL 
The following nonlinear system describes the dynamic model of the synchronous motor in the d-q 


reference frame [16], [17]: 
{ i(t) = f(x(t)) + g(x(t))u(t) + V(x) (1) 


where: 
w(t) —Fw(t) + SPA i (t) 0 0 
x(t) = fiat] ted = -BeO - Eih- o = [| esl 0 
ia(t) pu(t)ig — Bial) i 0 -ig 
a 
v(t) = | 0 
0 
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In which w is the rotor speed, (iq, iq) are the current components in the q — d axis, (ug, ua) are the 
stator voltage components in the q — d axis, Ly and Lg are the stator inductors in the q — d axis, Tm is the load 
torque (the load torque is an exogenous disturbance). The machine parameters are: the moment of inertia of 
the rotor J, the number of poles pairs p, the friction coefficient relating to the rotor speed f, the stator winding 
resistance R, the flux linkage of the permanent magnets À. In our study, the smooth-air-gap of the PMSM 
systems are considered, i.e., Lg = Lg = L. 


3. T-S FUZZY MODEL OF PMSM 
In order to represent the motor’s nonlinear model as a T-S model with the measurable (speed) param- 
eter as a decision variable. We rewrite (1) in the nonlinear state-space form shown B: 


l i(t) = A(w(t))e(t) + Bult) + DTm(t) a 
y(t) = Ca(t) 
where: 
"E P 9 oo) [-4 
A(w(t)) = |-™ R -pw(t)|,B=|¢ 0/,D=|0],Cc=[1 0 0] 
0 pw(t) -# 0 + 0 


Note that: w < w(t) < w [18], the nonlinear system of (2) could be represented by a T-S model with r = 2! 
fuzzy If-Then rules are as follows: 

Rule 1: IF z(t) is Fy, THEN «(t) = Ai a(t) + Biu(t) + DıiTm(t). 

Rule 2: IF z(t) is Fio THEN i(t) = Aga(t) + Bou(t) + DoT, (t). 

Where z(t) = w(t) denotes the premise variable, F}, and F32 refer to the membership functions which can be 
described as (3): 


Fu(w(t)) = OE, malt) =1- Faw) 3) 


the matrices of the local models are given by: 


f 3pA $ 3pà 1 

e a a oe, T a 

Aj = pe L pw ,Ag = y T pw , Bı = B2 = L 0 ‚Dı = D2 = 0 
0 po -È 0 pw -È 0 z 0 


based on singleton fuzzifier, the rule of product-inference, and the centre of gravity defuzzifier, the rules 
(1 and 2) become (4): 


i(t) = 2 hilz) (Aix (t) + Byu(t) + DiTm(t)) (4) 


where h;(z(t)) = EAO for all t > 0, hi(z(t)) > 0 and Y? h;(z(t)) = 1. 


4. FUZZY TRACKING CONTROLLER DESIGN 

The aim of this study is to design a fuzzy controller able to control the system’s state x(t) to track 
a desired set of VDVs x(t) and rejecting the impact of disturbance on the motor. The requirements for the 
feedback tracking control are as follows: x(t) — za(t) — 0 as t — oo from y(t) = y(a(t)), the desired output 
state yg(t) = y(xa(t)) is deducted. The tracking error is defined by Z(t) = x(t) — xa(t) and its time derivative 


is expressed by (5): 


Z(t) = i(t) — talt) (5) 


from (4) and the concept X ;—; hi(z(t))(Ai(x(t) — xa(t)), (5) will be {6}: 


&(t) = YO hi(z(t)){AiE(t) + Biu(t) + DiTin(t)} — talt) (6) 


i=l 
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In {6}, we present the new variable satisfying the following relationship: 


>D hi(z(t))BiT(t) = > hi(z(t))Biu(t) + +> hi(z(t)) Ayva(t) — ialt) (7) 


where, a new controller 7 (t) will be generated using the PDC method [12]. Using (8) tracking error system 
becomes: 


-Enel #(t) + Bir(t) + DiTm(t) (8) 


The tracking control problem is addressed in the following ways by the PDC controllers: 
Rule 1: IF z(w(t)) is Fj, THEN 7(t) = -Kı Z(t). 

Rule 2: IF z(w(t)) is Fig THEN 7 (t) = —K2Z(t). 

The final output of the fuzzy controller is determined by the summation: 


-= 2. hi(z(t)) K: (t) (9) 


applying control law (9) to model (8) leads the expression of the closed loop system: 


= Datel 2(t))(Ai — B:K;)ž(t) (10) 


i=1 g=1 


letting Gi; = (A; — B;K;), in becomes: 


z(t) = Dodo MEOH) GE) (11) 


5. STABILITY ANALYSIS 

The objective of this study is to develop a fuzzy state feedback controller in {9}, for system (11), able 
to driving the state of the system to track the desired variables. The gains of the PDC control law are 
determined using the subsequent theorem [13]. Theorem the equilibrium of closed-loop continuous fuzzy sys- 
tem (11) is asymptotically stable if there exist a symmetric matrix P. > 0, a diagonal positive definite matrix 
D. and matrices Q;; with: Qi; = Qii? and Qu = ore for i Æ j such that: 


GEP,+P.Gi+ D-P,D? <0,i=1,...,7 (12) 
(Git Giyrp ! p Gat Civ) Qy S0,8< 5 <0 (13) 
Qu . è . Qir 
>0 (14) 
Din. 2 5 Oe 


for i,j = 1, ..., r, s.t. the pairs (i, j) such that h;(z(t))hj(z(t)) = 0 Vt. 

The conditions of the previous theorem are related to a problem in the form of LMIs, that can be solved in an 
efficient way by convex optimization tools. The change of variable X = P71, K; = M;P7' and the use of 
congruence in inequalities to (14h, the LMIs expressions are given by: 


JX = X" S0,5%, = Y£ 3%, = Y} , Mi 


iin ji? 
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XAT + AX — BiM; — M] B]? +Ya XD? 


D eae (15) 

XAT + AX + XAT + A;jX — BiM; — M7 BY +2Yij < 0,i < j =1,...,7 (16) 
Yu Yio ... Yir 
Yis Yoo ... Yə 

a . | >0 (17) 
Yir Yor see Yrr 


we note that the controller gains K; that guarantees the tracking control is the same that stabilizes the system. 
The controller performance depends on the best choice of the diagonal matrix De. But, there is no specific 
method for choosing this matrix. Then, we propose here a new algorithm that calculates matrix D, and the 
controller gains based on the desired system characteristics. To this end, let us define the desired rise time Rig 
and the desired overshoot O,,g. The diagonal matrix is set to be in the following form: 


de 0 0 
0 -. 0) (18) 
0 0 de 


the proposed algorithm is summarized in the following diagram in Figure 2. 


Initial system parameters and D,(0) 


>. 
Vv 


Calculate the controller gains K; 
by resolving the LMIs (15), (16) and (17) 


Y 


Stabilization system using PDC approach 


de(i +1) = de (i) + step y 
7% calculate tr ,O,% of the system 


Figure 2. Proposed algorithm to calculate De and K; 


6. VDVS AND NONLINEAR CONTROLLER DESIGN 
As a way to identify the VDVs xq(t) and control law u(t) for the PMSM, (8) can be expressed as (19). 


2 hi(z(t))Bi(u(t) — 7)) = — 2 hi(z(t)) Azalt) — > hi(2(t))DiTm(t) + talt) (19) 


Assuming that: 


r r r 


g(x) = So hi(2(t)) Bi, A(a) = Y hi(2(t)) Ai pla) = So hi 2(O)Di (20) 


i=l i=1 i=l 
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in (19) can therefore be expressed in the compact form shown (21). 


g(x)((u(t) — r(t))) = A(@)@a(t) — 9(@)Tm(t) + talt) (21) 
Adapting the to the PMSM model, following matrix form is obtained (22). 


0 0 ig £ 3pA 0 Wd = i Wa 

1 Ug(t) — nl _ J Z, 7 a | 

T 0 P — Talt) = By L pw} |taa| = 0 | Tm(t) + A iqa (22) 
Lae 0 pw — $ ldd 0 idd 


Ta = [wa iqa iqaļ” is the vector of the VDVs. From (22), it remains that: 


B 3pA 1 
walt) = -F walt) + Fy galt) + 7Tm(t) (23) 
which induced that: 
2J. 1 
ialt) = -FE (Eal) + ealt) + 5Tm(t) (24) 


Noting that a flow model is not necessary for a PMSM. As a consequence, the angle of reference 


is determined by the position of the rotor. Additionally, because our machine has a smooth poles, the ideal 
value for its operation is determined when the internal angle is equal to 5, which indicates iaa = 0. Asa 


consequence, we were able to acquire the appropriate state vector shown (25): 


Wd B Ya 
talya, Tm) = |tga| = | Wat Gya + 5 Tm) es (25) 
tad 0 


where yq represent the required speed. From the second and the third expressions of (22), the nonlinear tracking 
control input is obtained: 
{ ug(t) = pAwa(t) + Riga(t) + LE + 14(t) (26) 


ualt) = —pLw(t)igact) + Talt) 


noting that the reference model obtained in depends on the desired speed and the load torque; this last 
latter is usually an inaccessible stat or immeasurable, that imposes to look for its estimated. In the following 
sections, we will show how to get the estimated load torque using both the extended Kalman filter and the 
disturbance estimator. 


7. EXTENTED KALMAN FILTER 

R.E. Kalman developed the Kalman filter in 1960. Due to developments in digital computer technol- 
ogy, the Kalman filter is a considerable study of research and application. Kalman filtering has been used in 
the field of manufacturing [19], navigation [20], aerospace [21], and others [22]-[24]. The EKF is an optimal 
recursive estimate approach for estimating the states of dynamic nonlinear systems based on the least-squares 
principle. To determine the conditional mean and covariance of the probability distribution of the state of a 
nonlinear stochastic system, it is therefore an ideal estimator by measurement noise and uncorrelated Gaus- 
sian process. The EKF may be used to estimate state variables since the state models are nonlinear. The 
back-electric and magnetic field (EMF) is assumed to be a state variable in this way. 


7.1. Problem formulation 
The following expression is a nonlinear discrete model using white noise: 


ee pr em (27) 
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where w and v represent the random disturbances. Knowing that v indicates measurement noise, which stands 
for measurement and sample errors and w refer to the process noise, which represents the parameter errors. 
The following definitions apply to the noise covariance matrices: 


Q =cov(w) = E {ww"} ,R=cov(v) = E {uv } (28) 


for the model’s linearization process, the partial derivative is adopted to give the discrete state models: 


| 
F(k) = | ——>—_ (29) 
o ( ðxT (k) z(k)=ê(k/k) 
error covariance matrix estimation: 
P~(k+1) = F(k)P(k)FT(k)+Q (30) 


where P~1(k + 1) represent a matrix of a priori error covariance. Calculating the gain of a Kalman filter: 
K(k +1) = P7 (k +1)HT(k)[H(k)P7(k+1)HT (k) + R] (31) 
the state estimation is given by: 
&(k +1) = &(k)K(k4+ 1)[y(k +1) — k(ê(k + 1))] (32) 
updating a matrix of error covariance: 


P(k +1) =(I1— K(k +1)H(k))P7(k +1) (33) 


8. EKF FOR THE PMSM 
The implementation of EKF for the PMSM requires the discrete motor model, which can be expressed 
as [19]: 


f 
a W—OO™m” 


rd. 
Te Pàis + (1-T;$)w Te Tm 
x(k + 1) = Tspwia JA (1 TÉ Jiq Tepi: } Ts Ly Uq (34) 
(1 — T, Ë )ia + Tspwiq + Ts -ua 


define matrix F and H as: 


1-T,4 T, “PA 0 

—Tsp(ia(k) + È) 1-Tsz#2 —Tspw(k) (35) 
Tspig(k) Tspw(k) 1-T.# 
w(k) 100 

h= |ig(k) = SM) = 0 1 0 (36) 
ig(k) stk) lo 0 1 


the EKF algorithm for the PMSM is given in Figure 3. 
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Initial Values 
Q, R, žo, Po 
] i 
( Prediction igi 
E(k +1) = f(&(k), u(k), Tm (k)) 
V P (k+1) = F(k)P(k)F("k) +Q Th, 
( Update Measurements N 
K(k +1) = P~ (k+ 1)H7 (k)[H(k)P—(k + 1)H" (k) + R] 
P(k +1) = (I — K(k + 1)H(k))P (k+ 1) —— 
&(k+1) =&(k) + K(k + 1)[y(k +1) — h(ž(k + 1))] Y 
P(k+1) =I- K(k +1)H(k)P (k+ 1) 


Figure 3. Extended Kalman filter algorithm 


9. LOAD TORQUE ESTIMATOR 

The reference model obtained in and the EKF algorithm shown in Figure 3 need the load torque 
state. However, this latter is generally an inaccessible state, or its measure is a difficult task that imposes 
looking for its estimated value using an approach based on the machine model combined with a proportional- 
integral controller as illustrated in Figure 4 [25]. The role of this controller is to eliminate the speed tracking 
error and guarantee the estimated load torque’s convergence to the load torque applied to the motor. The 
PI controller’s parameters can be determined using the pole placement technique. To this end, we write the 
mechanical equation of the PMSM as (37). 


dl eA O > TEE rs > K+ - sie 
Figure 4. Load torque estimator 
dw 
Te(t) — Tm (t) = J + fol) (37) 
Where T(t) is the electromagnetic torque which can be given by: 
3pr . 
T.(t) = ay iat) (38) 
which induced that: 
A A da 
T.(t) — Tn(t) = IF + FOE) (39) 
according to the scheme shown in Figure 4, the combined with a PI controller can be represented by the 
following matrix form: 
R (f+kp) kr D 3pà kp I 
W — tt -_-+ Ww See 2P 
a= J TI 233 J 4 
al-Pr E+E All w 
where x, (t) = [(@(t) — w(t))dt. 
The estimated load torque can be derived as follows: 
Tn = [kp bi] [2] +o -el (41) 
m p Le Pl day 
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by identifying with a second-order system, we obtain the following PI parameters: 
kp = 26,7 = f, ky = w2 J (42) 


€ and wpn denote the second-order system’s of damping ratio and the undimmed natural frequency. 


10. © SIMULATION RESULTS 

In this part, PMSM simulation tests have been conducted to validate the efficacy of the suggested 
technique. The overall control structure of the PMSM is shown in Figure 5, where the goal is to impose the the 
PMSM output speed to follow the desired speed. yg = wreg = 50rad /sec with a load torque Tm = 8N.m 
applied at t = 0.5s. 


SVPMW 


inventer 


iqa wa Ww 


T-S Fuzzy 
controller 


Figure 5. Overall control structure of the PMSM 


In order to facilitate the control data and clarify the results, an interface has been designed using 
Graphic User Interface (GUI); the interface has three figures: The first one is used to enter the different param- 
eters of the motor, the second one is used to determine the best choice of the diagonal matrix De. The third one 
is used to illustrate the simulation results according to the global diagram of the control structure diagram given 
in Figure 5. The GUI figures are shown in Figure 6 to Figure 8. The obtained controller gains that guarantee 
overshoot less than 1074% and time rise less than t = 0.005s are given by: 


K1i= 6.9313 16.2043 ee 


11.2537 14.0795 —0.1602 
0.6946 —0.2203 15.6902 


0.0886 —0.1153 12.9561 


Parameter values of the PMSM 


Stator resistance (Ohm) 1.67 


Longitudinal inductance (mH) 1.45 


Transverse inductance (mH) 1.45 


Stator Moment of inertia (Kg.rm’2) 0.013 


Number of poles pair 2 
| Friction coefficient (N.m.s/rad) 0.013 
Flux linkage (V. s/rad) 0.17 


Speed bound (rad/sec) -100| <we 100 


Validate Reset} 


Figure 6. The power curve of 5 MW wind turbine Figure 7. Performance coefficient as a function of 
Tip Speed Ratio 
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o5 0-5 4 
Simulation 
Simulate 
Clear 


L Stop ] 


xO ololo 6.9413 16.2043 -0.2969 


K1 0.6946 -0.2203 15.6902 


System initiale values 4 p Controller gains 
a olofo 


11.2537 140795 -01602 
EKF initiale Values——__ 2nd order system. K2 0.0886 -0.1153 12.9561 
Po [0.1 | 0.1 | 100 | 


Q |100 | 100 [1000 Das nia PI controller parameters: 


o6 
R Los | 0.1 | 0.3 5.4 


Figure 8. GUI figure to show simulation results 


The initial matrices and the EKF are chosen as follows: 


0.1 0 0 100 0 0 01 0 0 
W=] 0 01 0)],Q=] 0 100 O0 R=; 0 01 0 
0 0 100 0 O0 1000 0 0 O01 


The parameters of the PI controller are chosen as follows: Kp = 5.4,Kzr = 0.0202. The simulation re- 
sults of reference speed, actual speed and its estimates, reference q-axes current, actual q-axes current and its 


estimated, reference d-axes current, d-axes actual current and its estimated and the q-d control voltages are 
illustrated, respectively, in Figures 9 to 12. 


Speed (rad/sec) 
Current (A) 


10 0.5 0.505 0.51 0.515 0.52 4 
o 
o 0.2 0.4 0.6 0.8 4 o 0.2 0.4 0.6 0.8 i 
Time(sec) Time (sec) 


Figure 9. Speed reference, actual rotor speed and its Figure 10. Reference, actual and estimate currents of 
estimate q-axis 


Control voltages (V) 


0 0.2 0.4 0.6 0.8 1 


Time(sec) Time(sec) 


Figure 11. Reference, actual and estimate currents of 


Figure 12. Control input voltages 
d-axis 
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The simulation results demonstrate that the system is resilient to the disturbance with just a minor 
deviation during the time the disturbance is applied, the tracking error is obvious but very tiny, and the temporal 
response of the tracking is extremely low. The simulation results of load torque and its estimated, and the 
electrometric torque and its estimated are shown in Figure 13 and Figure 14, which demonstrate the good 
tracking. 


Torque (N.m) 


ie) 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 13. Actual load torque and its estimate 


10 T T T T 
eaei 


Torque(N.m) 


-2 


ie) 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 14. Actual electromagnetic torque and its estimate 


11. CONCLUSION 

This paper describes a novel fuzzy tracking control strategy for PMSM based on the T-S fuzzy model. 
Due to this, a fuzzy tracking control based on VDVs and the estimation of the unknown disturbance using an 
EKF and a DO has been proposed to develop a controller able to reject a completely unknown disturbance. 
Using Lyapunov’s approach, sufficient requirements for stability are determined. VDVs concepts have been 
utilized to simplify the construction of the reference model and control laws. A PMSM controller has achieved 
successful testing. The outcomes demonstrate that the suggested fuzzy tracking control approach may deliver 
the controlled system’s intended performances. 
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